Eur. Phys. J. C manuscript No. 

(will be inserted by the editor) 



Interacting dark fluid in the universe bounded by event horizon;A 
non-equilibrium prescription 



Subenoy Chakraborty Atreyee Biswas 



(N 

o 

ov 

(N 



43 

o 



Received: date / Accepted: date 

Abstract A non-equilibrium thermodynamic analysis has 
been done for the interacting dark fluid in the universe bounded 
by the event horizon. From observational evidences it is as- 
sumed that at present the matter in the universe is dominated 
by two dark sectors-dark matter and dark energy.The mu- 
tual interaction among them results in spontaneous heat flow 
between the horizon and the fluid system and the thermal 
equilibrium will no longer hold. In the present work,the dark 
matter is chosen in the form of dust while the dark energy 
is chosen as a perfect fluid with constant equation of state 
in one case and holographic dark energy model is chosen 
in the other. Finally, validity of the generalized second law of 
thermodynamics has been examined in both cases. 
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1 Introduction 

At present,based on the recent observational evidences 
(particularly from Type la supernovae observations[l])it is 
commonly believed that the matter in the universe is domi- 
nated by two dark components-dark matter( about 23%)and 
dark energy(above 73%). Dark matter(DM),the invisible mat- 
ter without pressure can explain galatic curves and large- 
scale structure formation while dark energy(DE),an exotic 
matter 
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with large negative pressure is responsible for the present ac- 
celerating phase of universe. Although there are several pro- 
posals for DE candidate,but still the nature of dark energy is 
completely unknown. 

However,the most common choice for DE candidate namely 
the cosmological constant A (having equation of state pa- 
rameter C0j\ = — l)has been discarded both from theoretical 
as well as from observational point of view. There is a huge 
order of discrepancy in the observed value and the theoret- 
ically estimated value of A [2] -there are two well known 
difficulties namely the "fine tuning" and the "cosmic coinci- 
dence" problems and are commonly known as cosmological 
constant problems [3]. 

However,there are different candidate for the dynamical DE 
scenario in the literature to interpret the present accelerat- 
ing phase of the universe namely a) the quintessence scalar 
field models[4] ,the phantom field[5],K-essence [6],tachyon 
field[7],quintom [8]etc. b)the DE models including Chaply- 
gin gas[9],brane world models [10], holographic and agegra- 
-phic[12] models and so on. 

From cosmological view point,it is interesting to consider 
interactions among the constituent matter components of the 
universe. But local gravity experiments put strong constraint 
on the interaction of DE with the baryonic matter[13] while 
there are no restrictions on the interaction among DE and 
DM;rather it is physically reasonable since DE gravitates - it 
may be accreted by massive compact objects(like BH,neotron 
star). But this energy flow from DE to DM should be small(but 
non-zero)from cosmological context. 

Initially,the coupling between DE and DM was considered 
to reduce the huge difference between theoretically predicted 
value and the observed value[14] of the cosmological con- 
stant and to solve the coincidence problem[ 1 5] .Further it has 
been shown that a proper choice of the interaction term may 
influence the perturbation dynamics and affect the lowest 
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multipoles of the CMB spectrum[16, 17]. Also recently,the 
analysis of the supernova data together with CMB and large- 
scale structure [18] revealed such interaction from expan- 
sion history of the universe. Further,in the context of the dy- 
namics of the galaxy clusters, signatures of the interaction 
between DE and DM has been analyzed[19]. 
Moreover,from thermodynamical view point the coupling 
between DE and DM has been studied [20] 
considering DE as perfect fluid with a well-defined temper- 
ature and it has been shown that at present epoch the energy 
flow should be from DE to DM for the validity of the second 
law of thermodynamics[21].In the present work, we consider 
the universe containing interacting DE and DM as the mat- 
ter constituent and it is assumed that universe bounded by 
the future event horizon is an isolated thermodynamical sys- 
tem. We shall consider thermodynamics of irreversible pro- 
cess (non equilibrium)due to energy flow and as a result the 
extensive property of the entropy of the whole system will 
no longer hold. We shall formulate the modified entropy of 
the whole system and examine the validity of the general- 
ized second law of thermodynamics. The paper is organized 
as follows: section 2 deals with a general prescription for ir- 
reversible thermodynamics with DE as a perfect fluid with 
constant equation of state while holographic DE model has 
been used in section 3. At the end in section 4,there is a brief 
discussion and concluding remarks. 



2 A study of the energy transfer between the dark 
sectors of the matter distribution: A general 
thermodynamic prescription 

In the context of the recent observational evidences,the present 
day accelerating universe is dominated by an interacting two 
fluid system-the dark energy(DE) and dark matter(DM).The 
energy conservation relations for both the subsystems are 



p m + 3Hp m = Q 
and 

p d + 3H(p d + Pd ) = -Q 



(1) 



(2) 



where the DE component is a perfect fluid having energy 
density and thermodynamic pressure p d and p d respectively 
while DM is in the form of dust having energy density p,„.The 
interaction term Q > indicates an energy flow from DE to 
DM. The explicit form of Q is chosen in the form [22] 



Q = 3HXp d 



(3) 



with A, a small dimensionless positive quantity. Although this 
choice of Q seems to be phenomenological but it is compat- 
ible with observations [22, 23] like SNIa,CMB, large scale 
structure,H(z),age constraints and recently in galaxy clus- 
ters. Also it has been shown to be effective in alleviating the 
coincidence problem. 



As we are considering the fluid system in the universe,com- 
-posed of two subsystems(DE+DM)at different temperatures 
interacting through exchange of energy,so it is reasonable to 
employ thermodynamics of irreversible process. 
Accordingly,starting from the Euler's relation: nTs =p+p, 
(n=number density of particles in a comoving volume and 
s= the entropy per particle) and using the above conservation 
relations (1) and (2) and the conservation relation 2 = —3H 
we have the evolution equations for temperature as 



^ = 3H- 
T m r 

and 



T d 



-^ = -3H{X + (o d ) 



(4) 



(5) 



where T,„ and T d are the temperature of the DM component 
and the DE subsystem respectively, (0 d (—l < (0 d < — I) is 
the equation of state parameter for the DE and r(= is 
the ratio of the energy densities of the two subsystems. Thus 
on integration we have 



T 



ro 



-{2+3(l+a> d )} 



and 



Td = M- 

v«0 



(6) 



(7) 



where Tq is the common temperature of the two subsystems 
in equilibrium configuration while ao,ro are the values of 
the scale factor and the ratio of the energy densities in the 
equilibrium state. It is to be noted that in deriving equation 
(6) one has to take into account of the temperature T,„ Q a~ 2 
for the DM sector in the absence of interaction. At very early 
stages of the evolution of the universe we have T m > T d and 
with the expansion of the universe,both the subsystems ap- 
proach to the equilibrium configuration with common tem- 
perature 7o(when a = ao)-Subsequently (i.e a > «o),the ther- 
mal equilibrium is violated due to a continuous transfer of 
energy from DE to DM with T m < Tq < T d .As we are con- 
sidering the universe bounded by the event horizon as an 
isolated system,so at the thermal equilibrium the common 
temperature Tq is nothing but the Hawking temperature at 
the horizon i.e. 7o = jEr^ wnere Re is the radius of the event 
horizon for the FRW model. 

For the present isolated system if we denote the entropies of 
the two subsystems as S m and S d and Se is the entropy of the 
bounding event horizon,then 



dS„ 



dQ„ 



dF 
dt 



dt dt 
and 

dS d dQ d dE d 



T d 



dt 



dt 



dt 



Pd 



dV 
dt 



(8) 



(9) 



Interacting dark fluid in the universe bounded by event horizon;A non-equilibrium prescription 



3 



while from the Bekenstein area formula, 
dS E 

-±=2kR e Re (10) 
at 

Here V = \ kR e is the volume of the universe bounded by 
the event horizon and E m = p,„V and E d = p d V. 
As the overall system is isolated so the heat flow across the 
horizon (Q/,) will satisfy 

Qh = -(Qm + Qd) (11) 
In equilibrium configuration,the entropy of the whole sys- 
tem depends on the energy densities and volume only and 
from the extensive property,it is just the sum of the entropies 
i.e. S,„ + S d + S E . However in non-equilibrium thermodynam- 
ics one has to take into account of the irreversible fluxes such 
as energy transfers in the total entropy and hence the time 
variation of the total entropy is given by 

dSr dS m dS d dSp ■ ■■ ■ ■■ 

-^ = ^ + ^ + ^-A d Q d Q d -A h Q hQh (12) 

where A d and A/, are the energy transfer constants between 
DE and DM within the universe and between the universe 
and the horizon respectively.Now using equations (8)-(10) 
the explicit form of different terms on the r.h.s of equation 
( 1 2) are given by 

\ 2+3(1+%) 



dS,„ 
dt 



dSj 
dt 

dS E 

~dT 



-3kH 2 r\ 



1 



-z J 



{1 



-3tzH 2 r\ 



£2 k -£2 d {l +XHR e )} 

3{X+a d ) 



1 



l+z 
co d + XHR E ) 



Q d x 



(1 

2kR e (HR e 



1 



AdQdQd 



A h QhQh 



-A d H«R E (l 

-HR E {3(1 
-H 2 R 2 E {3X 2 -\ 

\A d H*Rl{\ 



(O d ) 



co d 

2 



hXHR E )[2(l 
(3A-2)(1 + 



X(3(O d + q + l)}] 
+ Q k +CO d Q d ) [2(1 



f CO d ) 
co d ) + 3X] 



co d Q d HR E {3(l+co d ) 



+(O d n d )+2qHR E - 
+3X-2}] 
with z,the usual red-shift parameter. 

As the expression for c ^j- is very lengthy,so to get an idea 
about its sign we make use of the observed or estimated val- 
ues of different parameters present in the above expressions 
at present epoch as follows: 

(O d = -\,X = \,z eq = 5.56 x \0\r eq = 1.09 x \0 5 ,Q d = 
0.72, Q k = 0.02- z = 0,q = -0.57. 
So we have 
dSj 



dt 



= R E 5.9 x 10 5 (HR E ) 2 {HR E - 1.25)+ 

0.2A d {HR E )\l.92 -HR E ) +0.28A h (HR E ) 2 x 
(1.43 -HR E ) +6.2&(HR E - 1)] 



where A d =A d H 2 and A/, =A h H 



So if we take A d ,A h > and 1.25Ra < Re < IA3Ra we 
see that > i.e generalized second law of thermody- 
namics (GSLT)holds on the event horizon for the present 
irreversible thermodynamic al system. 



3 Holographic Dark Energy Model 

A typical dark energy model which satisfies the holographic 
principle is known as holographic dark energy (HDE)model. 
According to this model using effective quantum field the- 
ory the energy density is given by [11] 

Pd - -gr 

K E 

where 'c' is a dimensionless parameter which may be esti- 
mated from observation [24] and the radius of the event hori- 
zon is chosen as the IR cut-off length to obtain correct equa- 
tion of state and the desired accelerating universe [25]. So 
one can write Re as 

Re ~ ^k,H 

where Q d = is the density parameter. 
At first for simplicity of calculations we use the dark sector 
as the non-interacting two subsystems namely the HDE and 
the DM. Then the density parameter evolves as [16] 



Q' d 



Q d (l-Q d )[ 1 



2Ja n 



(13) 



and the variable equation of state parameter for the HDE is 



1 2JTT d 



«W = -3- 



3c 



(14) 



where '/' stands for differentiation with respect to x=lna. 
In thermodynamics, starting from Euler's relation,the tem- 
perature of the HDE (a perfect fluid with variable equation 
of state)can be written as 

T d = T cm {\ + (O d )e- i ^^ 

which on integration using (13) gives 



(15) 



T d = T d0 -(l-^)(l- 
ao c 



^^(l + y^ 



2 

x(l + -« 
c 



^T d y-o 



(16) 



with T d Q an integration constant and ao the value of "a" 
when DE and DM are in thermal equilibrium. Also the tem- 
perature of the DM subsystem (behaving as dust)varies as 
the reciprocal of the square of the scale factor[[22]]i.e 



T,„oa 



-2 



or, T,„ 



T m [ — 

,«0 



(17) 
(18) 
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Then by virtue of the extensive property,the entropy of the 
whole system is just the sum of the entropies of the subsys- 
tems and the entropy of the horizon,i.e, 

dSj dS m dS d dSE 
dt dt dt dt 
Then as before using Gibbs' law to obtain the explicit form 
of the first two terms on the r.h.s of above equation and using 
Bekenstein entropy-area formula for the 3rd term we obtain 



dSj 
dt 



2nR E 



HRe-1 



2 R E p m {l+Z eq f 
TmO (1+Z) 2 



4R E p d (l+z) 

37do(l +Zeq) 



(i-Vflrf)a*(i + VA,)^ 



(i 



-) 



Now using T m o = T d Q = j^-, the Hawking temperature as- 
sociated with event horizon when DE and DM are in equi- 
librium,we obtain 



^=2kR e (x-lx 2 ) 



where x = HR E , I = 4-n 



3(1-04)0 



2 I 2£i d b 



, a 



l+Ze, 
1+Z 



c I 



ft = (l-V5rf)^(l + VHrf) F3 (l 

Hence for the validi 
Re < PRa with a,j3 



Hence for the validity of GSLT we have / < | and ccRa < 



l Tv /l-4/ 
21 



system does not satisfy the extensive property (as in equilib- 
rium case),rather it modifies as equation (12). 
The temperature of the two dark sectors are now given by 

(2 2jO~ d 

T c i = T c [q — — A 

\3 3c J «o 

and 

r a 

where r = Qr.rn = §^,z = JThi and T m <\ and T d n are inte- 

Pd PdO 



gration constants and I = exp 



i f 



(l_,2 )(1 + ^ ) _3^ 7 2 



However if we restrict ourselves to flat universe,then 



Qm 

Qd 

T d 



QdQd 



QhQh 

dS E 
dt 



-1 



{\-O d ){\+ Zeq )I 
3nc 3 (l+z)H- 1 



1 



ro* 



Q' d 9c 4 H 

ol 4 



Q' d 9c 4 H 
Ol 4 



2 J n. 



3c 



Ic 1 

o^~o~ d 

' 2 

3c 



O d - 



^1 

9c 2 



2fl| 
9c 2 



2kR e R e = 6.2%cH~ 



C 

o~ d 



ol 

3c 



ol 

3c 



Xc\JO d (X 



■k — 



k 2 c 2 



3c 



1 

o~ d 



X + 3 



10, 



So the time variation of total entropy reads as 



Now we shall generalize our model by considering in- 
teraction between HDE and DM. The form of the interac- 
tion term is chosen same as in the previous section. Then the 
evolution of the density parameter and the equation of state 
parameter for HDE are modified as 

^^Od 



o' d = o d {i-o d )[\ 

and 

(0 d = 



-3l0 2 



<k + - 



I, 2JO, 



3c 



Now integrating the energy conservation equations the ex- 
plicit form of the energy density components are obtained 

as: 



Pm — PmO 

and 



— exp 



Pd 



a 

PdO [ — 

\«0 



exp 



3k I ** 
ar 



dz 



(l-z 2 )(l + |)-3Az 2 

In non equilibrium extended thermodynamics due to irre- 
versible fluxes like energy transfers the entropy of the whole 



dS T _ 3nc 3 H- 1 r 1 (l+z) 
dt 



1 



A-f + 



2jn. 



roVOd 



l-O d 



<Od 



lOd 



\ — A 



9c d 



2jSl n 



kc(k 



9c 4 Q' el H- 1 



1 

~O d 



2^2 



k 2 c 



201 

9c 2 



9c 2 



9c 




+6.2%cH- 1 ( — 
0,i 



IO„ 



with A =AH 2 

From the above expression for 4$jL it is not possible to exam- 
ine the validity of GSLT,however we give a graphical repre- 
sentation of --fr- with the variation of O d and I in figure 1 .The 
figure tells us that with the observed dominance of dark en- 
ergy(i.e O d > 0.75)GSLT is always satisfied. 
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dt 



5x10' 




Fig. 1 Graphical representation of '-jf- in the unit of H 1 with the variation of £2j and I for A 
3.56 x 10 7 ,z = 



\,z eq = 5.56 x 10' ',r eq ■■ 



1.09 x 10 , A,/ 



4 Discussion and Concluding remarks: 



A study of non-equilibrium thermodynamics for the uni- 
verse bounded by the event horizon has been done with mat- 
ter content as interacting two-fluid system - the two dark 
components known at present as dark matter and dark en- 
ergy.As usual,the dark matter is chosen in the form of dust 
while in two separate sections the dark energy is chosen as 
perfect fluid with constant equation of state and holographic 
dark energy model respectively.Irreversible thermodynam- 
ics is applied to the isolated system (i.e. universe bounded 
by the event horizon)as the mutual interaction between the 
two dark fluid species results in a spontaneous heat flow be- 
tween the horizon and the fluid system.At early epoch of 
the evolution of the universe the temperature of the DM is 
larger than that of DE and both approaches the Hawking 
temperature of the horizon in course of expansion of the 
universe. However,subsequently this equilibrium configura- 
tion is destroyed due to a continuous transfer of energy from 
DE to DM, and hence the extensive property of the entropy 
of the whole system can not be applicable to the present sys- 
tem.Though the expression for the time variation of the total 
entropy of the system is very complicated but it is possible 
to find restrictions for the validity of the GSLT in case of 
perfect fluid model of DE with present observed data and 
of holographic dark energy model (without interaction). In 
both cases radius of the event horizon is restricted in a range 
for which both the bounds are proportional to the radius 



of the apparent horizon. On the other hand,in case of holo- 
graphic dark energy with interaction,even the temperature 
can only be evaluated in integral form and hence no explicit 
analytic form for total entropy variation is possible. So ana- 
lytically we do not have any conclusion regarding validity 
of GSLT,however,we have shown only graphically 
Finally we make a comparison with previous similar work 
done by Karami et al [26] for universe bounded by apparent 
horizon. In that work they have only considered DE as per- 
fect fluid with constant equation of state and have shown the 
validity of the GSLT with a restriction using the present ob- 
served data.In the present work,our derivation of non-equili- 
-brium temperature of the two dark sectors with variable 
equation of state is new and also use of event horizon for 
interacting dark fluid is initiated in the present work.For fu- 
ture work, we shall try with other dark energy models for 
non-equilibrium thermodynamics. 
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